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A scaling consideration of the Kondo lattices is performed with account of logarithmic Van Hove 
singularities (VHS) in the electron density of states. The scaling trajectories are presented for 
different magnetic phases. It is demonstrated that VHS lead to a considerable increase of the 
non-Fermi-liquid behavior region owing to softening of magnon branches during the renormalization 
process. Although the effective coupling constant remains moderate, renormalized magnetic moment 
and spin- fluctuation frequency can decrease by several orders of magnitude. A possible application 
, to /-systems and weak itinerant magnets is discussed. 
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Tn , Anomalous rare-earth and sctinide compounds are studied extensively starting from the middle of 1980s [l|, |2|. 

They include so-called Kondo lattices (with moderately enhanced electronic specific heat) and heavy-fermion systems 

demonstrating a huge linear specific heat. Main role in the physics of the Kondo lattices 0-13 belongs to the interplay 

^ ' of the on-site Kondo screening and intersite exchange interactions. Following to Doniach criterion [5|, it was believed in 

early works that the total suppression of either magnetic moments or the Kondo anomalies takes place. However, 

later experimental data and theoretical investigations made clear that the Kondo lattices as a rule demonstrate 

magnetic ordering or are close to this. This concept was consistently formulated and justified in a series of the papers 

'^ i [a-Q treating the mutual renormalization of two characteristic energy scales: the Kondo temperature Tk and spin- 

P5 ' fluctuation frequency aJ. A simple scaling consideration of this renormalization process in the s — f exchange model 

O , [9, 10] yields, depending on the values of bare parameters, both the "usual" states (a non-magnetic Kondo lattice or 

a magnetic state with weak Kondo corrections) and the peculiar magnetic Kondo-lattice state. In the latter regime, 

small variations of parameters result in strong changes of the ground-state moment. Thereby high sensitivity of the 

CO ' ground-state moment to external factors like pressure and doping by a small amount of impurities (a characteristic 

^ ' feature of heavy fermion magnets) is naturally explained. 

Qy During 1990s, a number of anomalous /-systems (Ua;Yi_2;Pd3, UPt3_2,Pda;, UCu5_a;Pda;, CeCug-ajAua;, 

^I Ua;Thi_2;Bei3,CeCu2Si2,CeNi2Ge2, CerNia etc.) demonstrating the non-Fermi-liquid (NFL) behavior have become 

a subject of great interest (see, e.g., the reviews [ll|, ll2[). These systems possess unusual logarithmic or power-law 
temperature dependences of electronic and magnetic properties. 

Various mechanisms were proposed to describe the NFL behavior [1J|, including two-channel Kondo scattering 
[Ij, [13 , "Griffiths singularities" in disordered magnets [1^ , strong spin fluctuations near a quantum magnetic phase 
transition [13, [li- 
lt is important that experimentally the NFL behavior (as well as heavy-fermion behavior) is typical for systems 
lying on the boundary of magnetic ordering and demonstrating strong spin fluctuations [J, [ll| . 

The NFL behavior close to the quantum phase transition was theoretically studied in a number of works |18l42l] | . 

5^ In particular, renormalization group investigations near the quantum phase transitions connected with the topology 

5^ , of the Fermi surface were performed [20|] . A scaling consideration of the Kondo lattices with account of singularities 

in the spin excitation spectral function (which are owing to Van Hove singularities in the magnon spectrum) yields 

the NFL behavior in an extremely narrow interval of bare parameters only [9l| . As demonstrated in Ref. t22| , when 

taking into account renormalization of spin-excitation damping, the region can become considerably broader. 

The systems under consideration demonstrate both local-moment and itinerant features. Moreover, lar ge l inear 
specific heat and NFL behavior is observed also in some d-systems including layered ruthenates Sr2Ru04 [23| and 
Sr3Ru207 24]. 

It is well known that magnetism of itinerant systems is intimately related to the presence of Van Hove singularities 
(VHS) near the Fermi level. Therefore, it is instructive to treat the Kondo effect in systems with a singular electron 
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spectrum. This is the aim of the present paper. 

It is evident that the Kondo effect in such systems has a number of pccuhar features. In particular, for the 
logarithmically divergent density of states 



p{E) = A\n 



D 

B\E\ 



(the energy is referred to the Fermi level, D is the half-bandwidth, the constants A and B are determined by the band 
spectrum) the Kondo singularities at Ep become double-logarithmic. Perturbation expansion yields a non-standard 
expression for the one-centre Kondo temperature. 



Tk (xDcxp -1/{AI) 



Nl/2 



(1) 



instead of the result for the smooth density of states, Tk oc _Dexp[l/2/p(0)], / being the s — f exchange parameter. 
The logarithmic divergence in p{E) is typical for the two-dimensional case (in particular, for the layered ruthenates). 
However, similar strong Van Hove singularities can occur also in some three-dimensional systems like Pd alloys and 
weak itinerant ferromagnets ZrZn2 and TiBr2 [25, 26]. 

In the present work we consider the Kondo problem and the NFL behavior with the singular electron density of 
states for the lattice of d{f )-spms where a competition with spin dynamics takes place. 

In Sect. 2 the renormalization group equations in the presence of VHS are presented. In Sect. 3 the scaling behavior 
for a paramagnet and for magnetic phases with account of spin-excitation damping is considered. In Sect. 4 we treat 
the scaling behavior for the magnetic phases with account of the incoherent contribution to spin spectral function. In 
Conclusions, various electron properties and general physical picture of magnetism are discussed. 



II. THE SCALING EQUATIONS IN THE PRESENCE OF VAN HOVE SINGULARITIES 

We use the s — d{f) exchange model of a Kondo lattice 

H = Y^ ikcLck<T -^Yl ^*'^"/^4aCi/3 + Yl -^qS-qSq 

ko" iotj3 q 



(2) 



where tk is the band energy, S; and Sq are spin-density operators and their Fourier transforms, Jq are the intersite 
exchange parameters, a are the Pauli matrices. 

The density of states corresponding to the spectrum tk is supposed to contain a Van Hove singularity near the 
Fermi level. In particular, for the square lattice with the spectrum 



tk — 2t(cos kx + cos ky) + 4t'(cos kx cos ky + 1) 



we have the density of states 
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where K{E) is the complete elliptic integral of the first kind, the bandwidth is determined by |i^ — 8t'| < 4|t|. For 
t' — we derive 



piE) 




so that, according to ([T]), 



Tk cc D exp 
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Note that the expression for the Kondo temperature in the parquet approximation has a different form [27 
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Tk oc D exp 
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However, the expression ([5]) agrees with the numerical Wilson renormalization group calculation for the square lattice 
[28| , unlike the result (|5|; the corresponding problems of the parquet approximation in the Hubbard model are 
discussed in the works |29l |. 

In Refs.[3,[l3jl22l: the interplay of the Kondo effect and intersite interactions was investigated by the renormalization 
group method. This starts from the second-order perturbation theory with the use of the equation-of-motion method 
(within the diagram technique for pseudofermions, such an approximation corresponds to the one- loop scaling). 

We apply the "poor man scaling" approach [30[ . This considers the dependence of effective (renormalized) model 
parameters on the cutoff parameter C which occurs at picking out the singular contributions from the Kondo correc- 
tions to the effective coupling and spin-fluctuation frequencies. Using the results of Refs.[3, [23 we can write down 
the system of scaling equations in the case of the Kondo lattice for various magnetic phases. 

In the calculations below we use the density of states for a square lattice with t' — both in the two-dimensional 
(2D) and three-dimensional (3D) cases as a phenomenological one, so that 

e{E) = e{-D)F{E), g{-D) ^ ^, F{E) = 1° ^ + 1 (7) 

We adopt the definition of the effective (renormalized) and bare s — / coupling constant 

ffe/(C) = -2e/e/(C), g = -21 g, p = g{-D) (8) 

where C — ^ —0 is a flow cutoff parameter. Other relevant variables are the characteristic spin-fluctuation energy 
ZJefiC) and magnetic moment Sef{C). 

To find the equation for Ief{C) we have to treat the electron self-energy. For a ferromagnet (the case of an 
antiferromagnet is considered in a similar way, see Ref.Q) the second-order Kondo contribution reads 

^^1{E)^±2I^sJ: "'--'^ . (9) 

~^ ii — tk-q ± ^q 

where n^ — f{tk) is the Fermi function. Then we have 

«e/ = [4',)(£;)-Eg(i?)]/(25) (10) 

Picking out in the sums the contribution of intermediate electron states near the Fermi level with C < tk+q < C + SC 
we obtain 

«e/(C) - 2pFiC)I^7j{-^)6C/C (11) 

where r]{x) is a scaling function which satisfies the condition 77(0) = 1 which guarantees the correct one-impurity 
limit. In the magnetically ordered phase, oJ is the magnon frequency cjq, which is averaged over the wavevectors 
q =2k where k runs over the Fermi surface (for simplicity we use a spherical Fermi surface). In the paramagnetic 
phase (the problem of localized moment screening) ZU is determined from the second moment of the spin spectral 
density. 

Now we treat the singular correction to We/ and the effective magnetic moment Sef ■ We have within the spin- wave 
picture 

6S = -Y^Sib^b^) (12) 

q 

The singular contribution to magnon occupation numbers occurs owing to the electron-magnon interaction. Calcula- 
tion for a ferromagnet from the corresponding magnon Green's function yields [9| 

^ [ty, - tk-q - Wqj 

We see that, when considering characteristics of localized-spin subsystem, the lowest-order Kondo corrections orig- 
inate from double integrals over both electron and hole states. Then we have to introduce two cutoff parameters Ce 
and Ch with Ce~\- Ch = C (C is the cutoff parameter for the electron-hole excitations), 5Ce = —SCh to obtain 

dSef{C)/S = 2p^I^F{C/2)F{~C/2U-^)6C/C (14) 



The renormalization of spin-wave frequency owing to magnon-magnon scattering is given by 

(JWq/Wq = -aqS{blfiq)/S (15) 

Further on we pass to the magnon frequency averaged over the Fermi surface. Then we have (oq — > a) 

SuJef [C)/W = aSSef {C)IS = 2ap^l'^F{C/2)F{-C/2)ri{-'^^)5C/C (16) 

The latter resuh holds for all magnetic phases with a = 1 — a for the paramagnetic (PM) phase, a = 1 — a' for the 
antiferromagnetic (AFM) phase, a = 2(1 — a") for the ferromagnetic (FM) phase. Here a,a' ,a" are some averages 
over the Fermi surface (see Ref.j9|]), a' = in the nearest-neighbor approximation. This approximation enables us 
to use a single renormalization parameter, rather than the whole function of q. For simplicity, we put in numerical 
calculations below a = \ (although the deviation 1 — a just determines critical exponents for physical properties, see 
Ref.[23| and Conclusions). 

The scaling picture (which determines the NFL behavior) is influenced by not only real, but also by imaginary 
part of the spin- fluctuation energy. The latter is even dominating in the paramagnetic phase (e.g., in the Heisenberg 
model a spin-diffusion picture can be adopted at high temperatures). In the magnetically ordered phases, the damping 
comes from paramagnon-like excitations. They can be taken into account starting from the magnon picture of the 
localized-spin excitation spectrum. In the s — f exchange model the damping is proportional to I^ and to the magnon 
frequency (for the details, see Ref.[22]). The dependence of the damping on the magnetic moment S (which is strongly 
renormalized) is crucial for the size of the NFL region. The calculation of the damping in the second order in / yields 
the contributions of the order of both I^S and I^ [3l|, |33 (formally, they correspond to the first and second order 
in the quasiclassical parameter l/2S'V^The corresponding problems from a semiphemenological point of view are 
discussed in Ref. |33| . Similar to Ref. j22l| . here we do not introduce the additional factor of S to obtain in terms of 
renormalized quantities 

7e/(C) = kF{Cl2)F{-C/2)glf{C)uJ,f{C), (17) 

the factor k being determined by the bandstructure and magnetic ordering; we put in numerical calculations k — 0.5. 
When taking into account spin-wave damping 7 we have 

/ ^e/(C) \ ,^^ ^e/(C) 7e/(g) 



\C\ ) 'V \c\ ' \C\ 

Replacing in the Kondo corrections g — >■ gef{C), uj — >■ ZUe/(C) we derive the set of scaling equations with account of 
the energy dependence of the electron density of states: 

dgefiC)/dC = F(C)A, (18) 

d\nuJefiC)/dC = -aF{C/2)F{-C/2)A/2, (19) 

d\n'Sef{C)/dC = -F{C/2)F{-C/2)A/2 (20) 



with 



A = A{C\uefiC),-f,fiC)) = -^^V ( 1^1 , -^^- 



Similar equations can be obtained for the general SU{N) Coqblin-SchriefFer model [9| where a/2 — > a/N (there are 
some peculiarities for FM case owing to asymmetry of spin- up and spin-down states). 

One can see that the renormalizations of the spin- fluctuation energy ZJex (C) and the damping are more strong than 
that of gefiC) owing to the factors of F{±C/2). We obtain from ([19]), ^ 



S^fiC) _ fuJ^fiQY^^ ^21) 



S 
However, the simple expression for p — const, 

uJefiC) = ZJeM-a[9efiC) - g]/2), (22) 




FIG. 1: The scaling functions ri{x) for a ferromagnet (solid line), antiferromagnet (long-dashed line) and paramagnet (short- 
dashed line) in the 2D (a) and 3D (b) cases 



does not hold for the logarithmic density of states; no simple relation with the quantity gef{C) = F{C)gef{C) is 
obtained either. Expanding in l/ln|_D/C|, we derive 



ZO^fiO c. SJexp (~^ y"^ ^[gefiC) - g,f{C')] ~ aln2[g,f{C) - g] 



(23) 



Now we treat the scaUng functions rj. In the paramagnetic case we use the spin-diffusion approximation (dissipative 
spin dynamics) to obtain (cf. Ref.f^) 



"'"'<c) 



1 + I?(k - k') VC2 



, ZJ = AVkp 



(24) 



tfc=t.,=0 



where V is the spin diffusion constant, the averages go over the Fermi surface. Integration yields 



Pj\/, \ _ j arctanx/x d — i 

"^ W = I I 1 [1 ^ (1 ^ x2)l/2]/(l + .^2)|l/2 ^ ^ 2 

In the FM and AFM phases for simple magnetic structures we have 

7y(a7,//|C|,7e//|C|)=Re((l-(c.k_k,+z7k-k')VC^)"') 
For an isotropic 3D ferromagnet integration in (|25p for 7 = const and quadratic spin-wave spectrum Wq ex q^ yields 

(26) 



(25) 



■n 



--^.,z) = Aln(l + ^)' + ^' 



Ax (1 - a;)2 -I- z^ 



where x = ZUe//|C|, z — 7ef/|C|. Although details of the spin- wave spectrum are not reproduced in such an approach, 
the renormalization of spin-wave frequency (which is an cutoff for the Kondo divergences) is adequately reproduced. 
In the 2D case we obtain in the same approximation 



V^^^ix, z) = iRc{[(l + iz){\ + iz- x)]-i/2 + [(1 _ j^)(i _ i^ ^ x)]-^'^} 
For an antiferromagnet integration in (1251) with the linear spin-wave spectrum Wq ex cq gives 



V 



^^'Hx,z) 



2x 



^ {ln[(l + z^+ x^f - 4x2] „ 2 ln(l + z^)} d^3 
2Im{x^ - 1 - 2iz + z^)-^/^ d = 2 



(27) 



(28) 



which modifies somewhat the results of Ref. [9| . 

The plot of the functions r]{x) = t]{x,0) for different magnetic phases is shown in Fig.l. Note that for a 2D 
antiferromagnet r]{x) vanishes discontinuosly at x > 1. However, a smooth non-zero contribution can occur for more 
realistic models of magnon spectrum. 
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FIG. 2: The scaling trajectories (?e/(C ~ ^^\D/C\)/g (a) and the corresponding dependences ln(5'/S'(^)) — ln{uJ/uJe /{£,)) (b) 
for 2D paramagnets with the logarithmic density of states (solid lines, g = 0.029 < QcQ = 0.030 > Qc) and constant density of 
states (dashed lines, g — 0.14 < gc,g = 0.145 > gc), a = 1, A = ln(D/ZZJ) = 5 

III. THE SCALING BEHAVIOR IN PARAMAGNETIC AND MAGNETIC PHASES 



We start the discussion of scaling behavior from the simple case of the Coqblin-SchriefFer model in the limit A^ — > oo. 
Then the renormalization of the magnon frequency is absent and the scaling behavior can be investigated analytically, 
similar to the ^j. We have 



l/ge/(C)-l/.g-G(C) 






-D 



C" 
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7) 



The equation (|29p can be used even for N — 2 provided that g is considerably smaller than the critical value gc 
effective coupling g^fiC) begins to deviate strongly from its one- impurity behavior 



(29) 



The 



l/gef{C)c^l/g^-\n'\D/C\ 



(30) 



at \C\ 
G{C = 



The boundary of the strong coupling region (the renormalized Kondo temperature) is determined by 



-n) 



-1/g. Of course, T^ means here only some characteristic energy scale extrapolated from high 



temperatures, and the detailed description of the ground state requires a more detailed consideration. In the PM, 
FM and 2D AFM phases spin dynamics suppresses T^. To leading order in ln(Z?/ZU) we have 



n 



(Ti 



-2N1/2 



with Tk given by ([5]). [However, owing to the minimum of the scaling function ((28| (Fig. lb), in the 3D AFM case 
spin dynamics at not large ZU results in and increase of T^.] 

Provided that the strong coupling regime does not occur, i.e. g is smaller than the critical value gc, gefiC — > 0) 
tends to a finite value g* . To leading order in \n{D/uJ) we have 



'^2 uJ 2 



(31) 



(which yields also a rough estimate oi gc for N — 2). An account of next-order terms results in an appreciable 
dependence on the type of magnetic ordering and space dimensionality. For PM, FM and 2D AFM phases the critical 
value gc is given by l/ffc = — G(0), and in the 3D AFM case gc is determined by the minimum of the function G{C). 
The numerical calculations for A^ = 2 were performed for A — ln(_D/ZU) = 5. The plots are presented in Figs. 2-5 
for both smooth and singular bare densities of states. We compare these cases in the same relative interval of the 
coupling constant \g — gd/gc- The scaling process for finite N in the former case is described in Ref. Q- It turns out 
that the qualitative picture near the critical value of the coupling constant is rather universal for the same interval 
Iff ~ 9c\/gc depending mainly on the scaling function and the damping, but not on the details of the bare electron 
density of states (even on its singularities). The shift of the Van Hove singularity below the Fermi level does not 
influence strongly the results too, although we cross the singularity during the scaling process (the singularity is in 
fact integrable). 
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FIG. 3: The scaling trajectories QefiO/g {^) ^^'^ the corresponding dependences ln(S'/S'(^)) = ln{uJ/liJef{^)) (b) for 2D 
ferromagnets with the logarithmic (solid lines, g = 0.052 < gc,g ~ 0.053 > gc) and constant density of states (dashed lines, 
g = 0.186 < gc,g = 0.189 > gc), k — 0.5, other parameters being the same as in Fig. 2 
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FIG. 4: The scaling trajectories ge/(C)/S' (s-) and the corresponding dependences ln(5'/S'(^)) = ln(Zj/ZUe/(C)) (b) for 2D 
antiferromagnets with the logarithmic (solid lines, g = 0.0361 < gc,g ~ 0.0365 > gc) and constant density of states (dashed 
lines, g = 0.153 < gc,g = 0.155 > gc) 



An important quantitative difference in the presence of VHS is that the renormalized coupling constant is consider- 
ably smaller. Moreover, relative renormalization of the coupling constant is also smaller (Fig. 2-4). This makes using 
perturbation theory and lowest-order scaling analysis more physically reliable than for the smooth density of states. 

At the same time, in the presence of VHS the renormalization of spin-fluctuation frequency (and of magnetic 
moment) becomes larger (a "soft mode" situation, which favors a NFL behavior). Formally, this is due to that the 
derivative 51nZUe/(C)/5C is proportional to F^{C/2), and dgef{C)/dC to F{C) only [see Eqs. (Il9l)-(l20l)]. Such a 
situation is similar to the scaling in the large-/ limit {I is the number of scattering channels for conduction electrons) 
where the effective coupling is not renormalized, (2Z -I- \)g^ /2 = g^— const, since a -^ a{2l -I- 1) Q- 

For a paramagnet with pure dissipative dynamics, the one-impurity behavior l/5e/(C — ln|£^/C'|) = l/g — C is 
changed at ^ ~ A by a NFL-like (smeared quasi-linear) region where 

Inp/uJefm = aln[5/^(e)] ^ iaA/2)^ (32) 

with A < 2/a. Such a behavior takes place both for g < gc and g > gc in a wide region of S,, i.e. up to rather low 
temperatures (Fig. 2b, for a discussion of physical properties see Conclusions). 

For magnetic phases, the singularities of the scaling function rj{x — >■ 1) play the crucial role. A rather distinct NFL 
behavior takes place in a more narrow region where the argument of the function -q is fixed at the singularity during 
the scahng process, so that 



uJcf{C)c^\C\, \np/uJcf{0]^i- 
Then for a smooth density of states we obtain from ([22]) 

ge/(0-5-2(^-A)/a 



(33) 



(34) 




FIG. 5: The scaling trajectories gef{(,)/g (a) and the corresponding dependences ln(S'/S(^)) = ln{ijj/uJefiO) (b) for 3D 
antiferromagnets with the logarithmic (solid lines, g = 0.072 < gc,g ~ 0.073 > (?c) and constant density of states (dashed lines, 
g = 0.228 <gc,g = 0.231 > gc) 



However, in the presence of VHS the relation ([25)1 between ZuefiC) and gef{C) is more complicated. 

In the case of a constant magnon damping considered in Ref. |9] , the region p3|) is not too narrow only provided that 
the bare coupling constant g is very close to the critical value gc for the magnetic instability {\g — gc\/gc ~ 10~*-^10~^). 
However, when taking into account the magnon damping renormalization, this region is considerably wider, although 
being smeared [22|, and the influence of the energy dependence p{E) becomes stronger. 

One can see that VHS lead to a considerable increase of the effective moment S* = S{^ —> oo). For the chosen 
deviation \g — gd/gc ~ 1% the moment renormalization S/S* can make up one-two orders of magnitude; even for 
Iff ~ 9c\/gc ^ 10 — 20%, the moment can decrease by several times. 

Since in 2D antiferromagnet ri{x > 1) = 0, a sharp transition to the saturation plateau occurs, unlike the FM case 
(cf. Figs.3 and 4). 

The situation in 3D systems with a logarithmic density of states can be considered in a similar way. Surprisingly, 
for 3D ferromagnets to high accuracy the scaling trajectories (except for a narrow critical region near gc) and the 
critical values gc turn out to be very close to 2D case (for both smooth density of states and with VHS). Therefore we 
do not show the corresponding plot. On the other hand, for 3D antiferromagnets VHS does not lead to suppression 
of magnon frequencies (Fig. 5). This is due to the influence of the minimum in the scaling function (Fig. lb). 



IV. THE SCALING BEHAVIOR WITH ACCOUNT OF INCOHERENT CONTRIBUTIONS 

The magnon approximation used in the previous Section is not quite valid since this underestimates the role of the 
damping. In fact, the spin spectral function should have an intermediate form between large-damping and spin-wave 
pictures, containing both coherent (magnon-like) and incoherent contributions A simple attempt to construct the 
corresponding scaling function as a linear combination was performed in Ref. [9|. In particular, in the case of a 
ferromagnet we have near the magnon pole 



((^^l^"q)). 



2SZ^ 



7 + ((5^i^"q))r°' 



(35) 



where the inverse residue at the pole is determined by 



l/2q = 1 



/9IIqH 



(36) 



nq(a;) being the polarization operator of the magnon Green's function. Besides that, there exists the singular contri- 
bution which comes from the incoherent (non-pole) part of the spin spectral density. Then we get 



dgefiC)/dC = F{C)A 
d\nUcfiC)/dC = ~aF{C/2)F{~C/2)A/2 

d{l/Z)/dC = d\nSef{C)/dC = -F{C/2)F{-C/2)A/2 



(37) 
(38) 
(39) 




FIG. 6: The scaling trajectories gef{£,)/g taking into account incoherent contributions (a) and the corresponding dependences 
ln(S/S(^)) — ln(Zj/Zje/(^)) (b) for 2D ferromagnets with the logarithmic (solid lines, g — 0.0330 < gc,g ~ 0.0335 > pc) and 
constant density of states (dashed lines, g — 0.143 < gc,g = 0.145 > gc), k — 0.5, a — 1,X = ln(D/Zj) = 5 



where 



A = [glfiC)/\C\][Z7^,,^iu:,jiC)/\C\) + (1 - Z)r^,„eo/.(^e/(C)/|C|)] 



with r]coh — T] ■ The choice of rjincoh is a more difficult problem; here we put simply rjincoh — V 
According to ((39|) we have 
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1 + In = 
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(40) 



(41) 



Consequently, the increase of magnetic moment owing to the Kondo screening leads to a considerable logarithmic 
suppression of magnon contributions to the spectral density. 

The role of the incoherent contribution becomes important only provided that S/ S{^) and Z deviate appreciably 
from unity. However, such a moment suppression just occurs at passing the region of singularity in rjcoh{x = 1), the 
further scaling process being determined by the incoherent contribution. 

The corresponding scaling trajectories are shown in Fig. 6. Now we have a two-stage renormalization. One can 
see that the well-linear "coherent" behavior region (which is rather narrow in Fig. 6) is changed by a PM-like "quasi- 
linear" behavior p2p with increasing ^. This crossover occurs when the function ijcoh reaches its maximum value at 
the singularity (cf. Ref.!?!). 

The "quasi-linear" behavior, although being somewhat smeared, is pronounced in a considerable region of S, even 
for not too small I5 — 5c|. The difference with the PM case is in that gef{£.) increases considerably at the first stage 
of renormalization owing to the singularity of the function rjcoh- 

A similar consideration can be performed for the antiferromagnetic phase (cf. Ref.Q). The account of incoherent 
contribution results in a smearing of the non-monotonous behavior of (?e/(C) in the 3D AFM case, so that at small 
\g — gd the maximum in the dependences ge/(0 ^''^d \n\uj/uJe /{£,)] vanishes completely. 



DISCUSSION AND CONCLUSIONS 



In the general problem of metallic magnetism, the peaks in the bare density of states (which are usually connected 
with VHS) near the Fermi level play a crucial role. Here we have investigated their influence starting from the Kondo 
lattice (s — d{f) exchange) model. 

For g ^- gc we obtain the magnetic state with small effective moment S* and a NFL-type behavior. The corre- 
sponding dependences S{T) = Sef{\C\ — > T) describe an analogue of the "temperature- induced magnetism" j34{. 
Such a picture is based on the many-electron renormalization (compensation) of localized magnetic moments and 
differs outwardly from the ordinary mechanism for weak itinerant ferromagnets with small S, which are assumed to 
correspond to the immediate vicinity of the Stoner instability. 

However, the physical difference is not radical. In fact, a continuous transition exists between the highly-correlated 
Kondo lattices and the "usual" itinerant-electron systems. In particular, one may view Pauli paramagnets as systems 
with high Tji of order of the Fermi energy; for enhanced Pauli paramagnets like Pd, Pt, UAI2, where the Curie- 
Weiss holds at high temperatures, one introduces instead of the Kondo temperature the so-called spin-fluctuation 
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temperature. A combined description of the Kondo lattice state and weak itinerant magnetism has been considered 



recently by Ohkawa [35|- Remember that the Kondo systems with VHS near the Fermi level under consideration just 
possess high values of Tk (see the Introduction). 

In this context, it would be instructive to describe weak itinerant magnets not from the "band" point of view, but 
from the side of local magnetic moments which are nearly compensated. Since a number of cerium NFL systems 
demonstrate itinerant-electron behavior 13] and it is customary now to treat UPta, CeSix and CeRh3B2 as weak 
itinerant magnets, the second approach appears already by far less natural than the first (see Refs.|7l.[37|). From the 
formal point of view, perturbation calculations in the Hubbard model, which describes itinerant-electron systems, are 
similar to those in the s — d{f) model, provided that one postulates the existence of local moments. Besides that, for 
two-dimensional itinerant systems with strong spin fluctuations the semiphenomenological spin-fermion model can be 
used which separates electron and spin degrees of freedom and is somewhat similar to s — d exchange model [36, 38:]. 

Further on, the question arises about the role which many-electron effects play in the "classical" weak itinerant 
3D magnets like ZrZn2 and TiBe2. Indeed, one can hardly believe that the extremal smallness of S in these systems 
is due to accidental bare values of N{Ef) and Stoner parameter. Moreover, the Stoner criterion is not valid even 
qualitatively (in particular, due to spin fluctuations the critical coupling Uc in the Hubbard model is finite when the 
Fermi level tends to VHS [38i]). Thus a scaling consideration in the presence of VHS would be of interest, especially 
with account of chemical potential renormalizations (cf. Ref.[23]). In particular, owing to VHS the chemical potential 
can depend weakly on electron concentration (the pinning phenomenon [39|). In the Kondo systems such a treatment 
may lead to a renormalization of the scale I5 — 5c| itself. 

Now we discuss some real layered systems. Specific heat is considerably enhanced in ruthenates Sr2Ru04. A 
gradual enhancement of the electronic specific heat and a more drastic increase of the static magnetic susceptibility 
were observed in Sr2-jyLayRu04 with increasing y. Furthermore, the quasi-2D Fermi-liquid behavior observed in pure 
Sr2Ru04 breaks down near the critical value y = 0.2. The enhancement of the density of states can be ascribed to 
the elevation of the Fermi energy toward a Van Hove singularity of the thermodynamically dominant Fermi-surface 
sheet. The NFL behavior is attributed to two-dimensional FM fluctuations with short-range correlations at VHS j23| . 

The bilayered ruthenate system Sr3Ru2 07 in the ground state is a paramagnetic Fermi liquid with strongly enhanced 
quasiparticle masses. The Fermi-liquid region of the phase diagram extends up to 10-15 K in zero field and is 
continuously suppressed towards zero temperature upon approaching the critical field of _B = 8T. In the vicinity of 
the putative quantum critical end point, NFL behavior has been observed in various macroscopic quantities including 
specific heat, resistivity and thermal expansion and has been described on the basis of phenomenological models 24] . 

We can mention also some layered /-systems. The layered Kondo lattice model was proposed for quantum critical 
beta-YbAlB4 where two-dimensional boron layers are Kondo coupled via interlayer Yb moments [40|. CeRuPO seems 
to be one of the rare examples of a ferromagnetic Kondo lattice where LSDA+U calculations evidence a quasi-2D 
electronic band structure, reflecting a strong covalent bonding within the CeO and RuP layers and a weak ioniclike 



bonding between the layers [41 1 



To describe layered antiferromagnetic cuprates, the 2D t — t' Hubbard model is often used which also describes 
Fermi-liquid and NFL regimes. Despite the density-of-state logarithmic singularity, the staggered spin susceptibility 
in this model does not diverge within the Fermi-liquid approach, the reason being the appearance of the logarithmic 
singularity in the quasiparticle mass [42|. The effective mass renormalization is beyond our lowest-order (one-loop) 
scaling consideration, but may play a role at an accurate treatment. The two-loop considerations of the flat-Fermi- 
surface and i — t' Hubbard models 43, 44] yield an (generally speaking, anisotropic) suppression of the quasiparticle 
weight (inverse effective mass) along the Fermi surface, the staggered spin susceptibility remaining divergent, although 
the divergence is considerable weakened. 

Note that the picture in the Hubbard model and s — d(f) model (where "direct" exchange interaction is present) 
can be considerably different. Recently, the e-expansion has been used for a scaling consideration of the 2D antifer- 
romagnetic Kondo lattice with the use of non- linear sigma model |45| . 

In the case of a smooth electron density of states, various physical properties of NFL systems are discussed within our 
approach in Ref. [22] . The temperature behavior of magnetic characteristics S and uj, which depend exponentially on 
the coupling constant, is decisive for the NFL picture under consideration. At the same time, the presence of VHS near 
the Fermi level influences strongly all the electronic, magnetic and transport properties in the scaling approach, as well 
as in the one-electron theory. The replacement p^g^JT) — >■ p^{T)glAT) with p{T) being considerably temperature 
dependent owing to VHS may modify somewhat the behavior of observable quantities. 

Consider the temperature dependence of the magnetic susceptibility x '^ S/uJ. Using the scaling arguments we can 
replace uJ — >■ a7e/(C), S — >■ Sef{C) with \C\ ~ T, which yields x(^) °^ T~'^ . The non-universal exponent Q is determined 
by details of magnetic structure (the difference a — 1 can be used as a perturbations, see Ref. [24I). Besides that, a 
number of crossovers are characteristic for NFL behavior under consideration. In the coherent "magnon" regime we 
have C, = [a— l)/a, and in the "quasilinear" (incoherent) region C = (a — 1)^/2. 

The temperature dependence of electronic specific heat can be estimated from the second-order perturbation theory, 
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Cei{T)/T (X \/Z{T) where Z{T) is the residue of the one-electron Green's function at the distance T from the Fermi 
level. Then we have 

C,i{T)/T (X p\T)glf{T)S,f{T)/uJ,f{T) « p\T)glf{T)x{T) (42) 

which results in a non-trivial behavior of the Wilson ratio. 

Following to Ref. |22| . a simple estimation of the transport relaxation rate (which determines the temperature 
dependence of the resistivity owing to scattering by spin fluctuations in AFM phase) yields 

- ex T^p\T)gl^{T)Sef{T)/uJe.{T) ex T^Cei{T)/T (43) 

However, a more refined treatment in spirit of Ref. [42] would be useful in some cases. 

In Ref. [22.], a mechanism of NFL behavior owing to peculiar behavior of spin spectral function was proposed. Here 
we treated a similar, but somewhat more simple and natural mechanism which is connected with the singularities in 
the bare electron spectrum. Of course, a more accurate treatment of magnetic fluctuations near the quantum phase 
transition is required. Therefore detailed investigations of the NFL behavior for a realistic Fermi surface and spin 
spectral function are of interest. An accurate investigation of the situation where VHS is shifted from Ep or two 
peaks are present below and above Ep |25|] would be also instructive. 

The research described was supported in part by the Program "Quantum Physics of Condensed Matter" from 
Presidium of Russian Academy of Sciences. The author is grateful to M.I. Katsnelson and A. A. Katanin for discussions 
of the problem. 
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